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THE  SOLUTIONS  OF  NON-HOMOGENEOUS  LINEAR  DIFFERENCE 


EQUATIONS  AND  THEIR  ASYMPTOTIC  FORM* 

BT 

K.  P.  WILLIAMS 

Introduction. 

Very  important  progress  has  recently  been  made  in  the  analytic  theory  of 
homogeneous  linear  difference  equations. 

Galbrun  f  has  used  the  Laplace  transformation  to  derive  important  ex- 
istence theorems,  and  has  investigated  the  nature  of  certain  principal  solu- 
tions for  large  values  of  the  variable. 

At  about  the  same  time  Norlund  J  applied  the  theory  of  factorial  series 
to  linear  difference  equations,  showed  the  existence  of  the  same  solutions, 
and  gave  their  asymptotic  form.  In  addition  to  the  case  of  polynomial  co- 
efficients he  has  considered  the  case  where  the  coefficients  can  be  expressed 
in  factorial  series.  In  a  more  recent  paper  he  has  investigated  an  equation 
with  still  more  general  coefficients.! 

By  means  of  a  method  of  successive  approximation,  and  a  suitable  extension 
of  a  contour  integral  due  to  Guichard,  Carmichael  ||  has  independently 
shown  the  existence  of  solutions,  and  has  found  the  bounds  of  increase  and 
decrease  of  the  solutions  in  a  direction  parallel  to  the  real  axis. 

The  latest  important  contribution  to  the  general  theory  is  by  Birkhoff.^I 
He  employs  a  matrix  notation,  and  shows  in  a  direct  manner  the  existence 
of  certain  intermediate  solutions  and  of  the  principal  solutions.  The  asymptotic 
form  of  these  solutions  is  determined  throughout  the  complex  plane.  A 
modification  of  the  integral  used  by  Carmichael  plays  an  important  role. 

*  Presented  to  the  Society  February  24,  1912. 

tComptes  Rendus,  vol.  148  (1909),  p.  905;  vol.  149  (1909),  p.  1046;  vol.  150  (1910), 
p.  206;  vol.  151  (1910),  p.  1114;  Acta  Mathematica,  vol.  36  (1912),  pp.  1-68. 

tComptesRendus,  vol.  149  (1909),  p.  841;  Bidrag  til  de  lineaere  Dijfernsligningers 
Theori.    Dissertation,  Copenhagen,  1910. 

§M6moires  de  l'Academie  Roy  ale  des  Sciences  et  des  Lettres 
de  Danemark,  Copenhagen;  7th  series,  vol.  6  (1911),  pp.  309-327. 

||  These  Transactions,  vol.  12  (1911),  pp.  99-134. 

«H  These  Transactions,  vol.  12  (1911),  pp.  243-284. 
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He  finds  also  the  form  of  the  periodic  functions  which  connect  the  principal 
solutions,  and  formulates  for  the  first  time  the  properties  which  are  sufficient 
to  characterize  a  given  system. 

While  the  analytic  theory  of  homogeneous  linear  difference  equations  has 
thus  been  extensively  treated,  no  general  theory  has  been  developed  for  non- 
homogeneous  equations,*  although  a  number  of  equations  of  particular  form 
have  been  considered  (see  Carmichael,  loc.  cit.,  Introduction). 

The  purpose  of  the  present  paper  is  to  prove  the  existence  of  certain  principal 
solutions  of  the  system  of  non-homogeneous  equations 

n 

gi  (  x  +  1 )  =  53  0>ij  (  »  )  9i  (  * )  (i  =  l,2,  —  .n-l), 

(1) 


9 


n{x+  1 )  =  X)  anj  ( X )  Qi  ( x )  +  b  ( X  +  1 ) , 


and  to  develop  their  asymptotic  form.  The  functions  a#  (x)  and  b  (x)  will 
be  assumed  to  be  rational,  f  with  poles  of  order  at  most  fx  and  v  respectively 
at  infinity,  so  that  we  can  write 


aij(x)  =  a^Otf  +  ^H  ) 

b(x)  =  x"(b+  ^+---) 


for  I  x  I  >  R , 


On  account  of  the  gain  in  simplicity  of  the  formulae  we  have  taken  a  system 
where  all  but  one  of  the  equations  are  homogeneous.  It  is,  furthermore,  evident 
that  the  system  (1)  would  be  no  more  general  if  all  the  equations  were  non- 
homogeneous,  since  the  solution  of  such  a  system  could  be  obtained  by  adding 
together  the  solutions  of  the  n  systems,  of  the  above  type,  which  would  be 
obtained  by  striking  out  the  non-homogeneous  term  from  all  the  equations 
except  the  first,  second,  .  .  .  ,  last,  in  turn.  Thus  the  theorems  for  (1)  can 
be  extended  readily  to  a  system  of  n  equations,  all  non-homogeneous. 

The  method  employed  is  based  on  a  method  analogous  to  that  of  varia- 
tion of  constants  in  the  theory  of  linear  differential  equations.  This  method 
has  previously  been  supposed  to  yield  only  formal  results. J    It  will  appear, 

*  It  is  possible  to  reduce  a  non-hornogeneous  equation  to  a  homogeneous  equation.  The 
application  of  the  general  results  for  a  homogeneous  equation  will  show  the  existence  of  solutions, 
but  gives  no  direct  means  of  studying  their  properties.  This  method  is  suggested  by  Wal- 
lenberg and  Guldbebg,  Theorie  der  linearen  Differenzengleichungen,  p.  87. 

f  Most  of  the  results  which  will  be  obtained  apply  with  slight  modification  to  the  case  where 
an  ( x )  and  b  ( x )  have  the  form  of  rational  functions  only  at  infinity. 

t  Wallenberg  und  Guldberg,  p.  vii,  "...  die  Methode  der  Variation  der  Konstanten, 
die  wegen  der  mit  ilir  verbundenen  Summationen  meist  nur  formale  Bedeutung  hat .  .  .  ." 
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however,  that  in  general,  if  not  in  all  cases,  the  sum  formulae  that  arise  can  be 
interpreted  so  as  to  give  solutions  analytic  in  the  finite  plane  except  for 
poles. 

There  is  a  difference  of  treatment  according  as  /x  >  0 ,  fj.  <  0 ,  or 
H  =  0. 

In  §  2  the  reduction  of  the  problem  to  that  of  making  valid  one  of  four 
formal  types  of  sums  is  effected. 

In  §  3,  under  the  hypothesis  that  ix  >  0 ,  a  direct  summation  to  the  right 
is  used  to  derive  a  first  principal  solution  of  (1),  which  we  denote  by  gn  (x) , 
#21  (x) ,  •  •  •  ,  <7wi  (x)  ■  This  solution  is  analytic  in  the  finite  plane  save  for 
poles  at  determinate  points  near  the  negative  real  axis.  In  §  4  it  is  proved 
that  in  the  right  half  plane  this  solution  is  asymptotic  to  the  series  formally 
satisfying  (1). 

It  is  shown  in  §  5  that  the  formal  solution  of  (1)  obtained  by  a  direct  sum- 
mation to  the  right,  such  as  is  used  in  §  3,  is  independent  of  the  fundamental 
set  of  solutions  of  the  homogeneous  system  associated  with  (1)  which  is  used 
to  determine  it.  This  principle  is  used  in  §  6  to  determine  the  asymptotic 
form  of  the  first  solution  in  the  left  half  plane. 

A  summation  by  means  of  a  contour  integral  to  the  left  is  employed  in  §  7 
to  determine  a  second  principal  solution,  gn  (x) ,  g2o  ( x ) ,  •  •  •  ,  gn2  ( x ) , 
analytic  save  for  poles  at  determinate  points  near  the  positive  real  axis.  This 
solution  is  asymptotic  to  the  formal  series  in  the  left  half  plane. 

In  §  8  the  relation  between  the  principal  solutions  is  determined. 

It  is  seen  in  §  9  that,  if  n  <  0 ,  the  same  theory  holds  except  that  the 
roles  of  the  right  and  left  half  planes  are  interchanged.  Similar  results  are 
proved  in  §  10  for  the  intermediate  case  /x  =  0 ,  except  in  certain  extreme 
cases. 

An  application  of  the  general  theory  is  made  in  §  11  to  a  single  equation 
of  the  nth  order. 

The  writer  wishes  to  express  his  appreciation  and  thanks  to  Professor 
Birkhoff  for  the  valuable  suggestions  and  aid  that  have  been  received  from 
him. 

§  1.  The  fundamental  existence  theorems  for  a  homogeneous  system. 
The  system 

(2)  gi(x+  1)  =  Zo(,(i)ft(i)  «  =  i,  2,  ....  n) 

is  the  homogeneous  system  associated  with  (1).  As  we  shall  make  use  of  its 
solutions,  it  will  be  convenient  to  enumerate  some  of  their  properties.  We 
state  below  the  existence  theorems  for  (2)  as  they  are  given  by  Birkhoff. 
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There  exist  in  general  n  sets  of  series: 

*y(iO  =  X'Me(^<T")'»r'(*W  +  ^+  ••'), 
(3)  s2j(x)  =  x»x(Pje-»yxr>ls2j  +  -^-\  J,    0'  =  1,  2,  n), 


snj(x)  =  x**(p,e-naxrj{sn,-  +  ^+-..), 

each  of  which  will  formally  satisfy  (2),  but  which  in  general  diverge. 
The  quantities  pj  are  the  roots  of  the  determinant  characteristic  equation 

(4)  |  a{j  —  8ij  p  |  =  0, 

where  5#  =  0  for  i  4=  j ,  and  8 a  =  1 ,  and  they  are  so  ordered  that 

(5)  I  Pi  [  5»  1  P2  I  ^  *  •  '  ^  I  Pn  |  • 

Furthermore,  they  are  all  supposed  different  from  zero,  that  is 

(6)  A—\  a{j  |  =  pi  p2  •  •  •  p„  H=  0 . 
It  is  assumed  that  such  series  exist,  and  also  that 

(7)  S=K;l*0. 
We  can  now  state  the  theorems  we  shall  use. 

There  are  two  fundamental  sets  of  solutions  of  (2),  the  first  and  second  principal 
solutions;*  they  may  be  denoted  by  hij  (x) ,       (x) ,  hnJ-  (x)  (j  =  1, 

2,  ft),  and  hj  (x) ,  h2}-  (x) ,  hnj  (x)  (j=l,  2,  ft),  re- 

spectively. They  have  the  defining  property  that  in  any  right  half  plane 
ha  (x)  ~  Sij  (a;),t  while  in  any  left  half  plane  h^  (x)  ~  s#  (x)  (i,  j  =  1 , 
2  ,  •  •  •  ,  n) .  The  functions  h^  ( x )  are  analytic  in  the  finite  plane  save  for  poles 
at  the  poles  of  ( x ) ,  and  the  zeros  of  the  determinant  \  a^  ( x )  | ,  or  points 
congruent%  to  these  points  on  the  left;  while  the  functions  hij  (x)  are  analytic 
save  for  poles  at  the  poles  of  aij  (x) ,  and  points  congruent  on  the  right. 

While  all  the  results  which  will  be  derived  could  be  obtained  by  making 
use  of  only  the  principal  solutions  of  (2),  it  will  simplify  the  work  to  make 

*  Our  first  and  second  principal  solutions  are  the  second  and  first  principal  solutions,  respec- 
tively, as  given  by  Birkhoff.    For  the  present  purpose  it  is  convenient  to  take  them  as  above. 

f  This  means  that  for  any  c  the  difference  between  ha  (x)  and  k  +  1  terms  of  s#  (  x  )  can 
be  written  x*x  {pj  e-v-  )xxrJ  $k  (x) /xk ,  where  lim  f4  (x)  =  0,  uniformly,  for  —  57r3=arg  (x  —  c) 

x=<x> 

=  \-k  or,  in  a  special  case,  for  —  \it  —  e  =  arg  (x  —  c )  =  fir  —  e ,  e  arbitrarily  small  (see 
Galbrdn,  loc.  cit.).  We  say  then  that  hi,(x)  is  asymptotic  to  Sij{x),  or  in  symbols 
ha  ( x )  ~     (  x  ) ,  and  call     ( x )  the  remainder  coefficient. 

t  The  points  •••  x  —  2  ,  x  —  l,x,x  +  l,x-{-2,  •••  are  called  a  congruent  set  of  points. 
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use  of  the  properties  of  one  of  the  two  sets  of  intermediate  solutions  that  are 
given  by  Birkhoff  (loc.  cit.  §  3). 

There  exists  a  set  of  intermediate  solutions  h\j  ( x ) ,  h'2j  ( x ) ,  •  •  •  ,  h'nj  ( x ) 
( j  =  1,2,  •  •  •  ,  n) ,  defined  only  at  a  sufficient  distance  from  the  axis  of  reals. 
The  elements  of  these  solutions  are  analytic  in  the  region  where  they  are  defined, 
and  are  such  that  h'a  (x)  ~  Sij  (x)  ivith  regard  to  x  in  the  left  half  plane*  and 
with  regard  to  v  (x  =  u  +  v  V  —  1 ) ,  in  any  right  half  plane,  f 

Although  these  intermediate  solutions  are  not  defined  near  the  real  axis, 
they  will  be  of  service  to  us  on  account  of  the  fact  that  they  maintain  their 
asymptotic  form  with  regard  to  v  in  any  right  half  plane,  and  because  the 
periodic  functions  which  connect  them  and  the  second  principal  solutions  are 
of  a  simple  nature.  The  relation  between  these  intermediate  and  the  second 
principal  solutions  is  given  by  the  formula 

hi,  (x)  =  cpij  (x)  h'n  ( x )  +  <p2j  (x)  h'i2  (x )  + 

(8) 

+        j-  (x)  hi,  i-i(x)  +  hi,  (x)   (i,  j  =  1,2,  •■•,7i), 

where  <p\,  (x) ,  cp2j  (x) ,  •  •  ■  ,  <pj-i,  j  (x)  are  periodic  of  period  1 ,  and  are 
analytic  at  a  sufficient  distance  from  the  real  axis.J 

When  we  put  j  =  1  in  (8)  we  see  that  the  solution  hn  (x) ,  •  ■  • ,  hnl  (x) 
(i.  e.,  the  one  corresponding  to  the  largest  |  p  \ )  maintains  its  asymptotic 
form  in  the  right  half  plane  with  regard  to  v .  These  functions  form  the  ele- 
ments of  the  first  column  of  the  determinant  \hij  (x)  \.  It  is  also  proved 
in  Birkhoff 's  paper  that  |  {x)\,  or  any  minor  formed  from  its  first  k 
columns  (k  <  n) ,  is  asymptotic  to  \sa  (x)  \,  or  its  corresponding  minor, 
with  regard  to  v  in  the  right  half  plane. 

There  exists  a  second  set  of  intermediate  solutions  h\j  ( x ) ,  such  that 
h'ij  ( x )  ~  (x)  with  regard  to  x  in  the  right  and  with  regard  to  v  in  the  left 
half  plane. 

§  2.    The  formal  solution  and  the  operation  of  summation. 

We  are  thus  assured  of  the  existence  of  solutions  of  a  homogeneous  system. 
This  will  enable  us  to  determine  easily  a  formal  solution  of  (1).  For  this 
purpose  let  hj  (x) ,  hj  (x) ,  •  •  ■  ,  lnj  (x)  (j  =  1 ,  2,  •  •  •  ,  n)  be  any  funda- 
mental system  of  solutions  of  (2).  Write 

n 

(9)  gi  ( x )  =  S ( a; )  la  ( x )         (i  =  l,  2,  ••♦,»), 

*  That  is,  asymptotic  in  the  sense  already  defined. 

f  This  means  that  the  difference  between  hij  (x)  and  k  +  1  terms  of  s;;-  ( x )  can  be  written 
x»x  ( pj  e-v  )x  xri  tk{x)  v~k  where  lim  £k(x)  =0,  —  iir  ^  arg  (x  —  c)= 

v=co 

J  Birkhoff,  loc.  cit.,  p.  266. 
Trans.  Am.  Math.  Soc.  14 
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and  let  us  determine  the  quantities  coi|(x) ,  co2  (x) ,  •  •  •  ,[co,;  (x)  so  that  <7i(x) , 
g2  ( x ) ,  •  •  •  ,  <7„  ( x )  will  be  a  solution  of  the  non-homogeneous  system.  If  we 
substitute  the  latter  expressions  in  (1),  after  first  putting 

COy  (  X  +  1  )  =  COy  (  X  )  +  AcOy  (  X  )  , 

we  find  in  virtue  of  (2) 

LAco,-  (x)  Uj  (x+  1)  =  0  (i  =  1,  2,[--.  ,  n-1), 

i=i 

EAcOy(x)Z„y(x+l)  =  6(X+  1). 
J'=i 

This  is  a  system  of  to  linear  non-homogeneous  equations  in  the  n  quantities 
Acoi(x),  Aco2(x),  Aco„(x).    Let  L  (x)  represent  the  determinant 

]  Uj  ( x )  | .  Further,  let  Lj  ( x )  be  the  minor  of  the  element  in  the  last  row 
and  jth  column  of  L  (x) .    From  (10)  we  have 

(11)  A«y  (  X  )  =  (  -  1  )«*»  £    . 

We  are  thus  naturally  led  to  the  consideration  of  the  simple  non-homogeneous 
equation 

(12)  A/(x)  =/(a+  1)  -  f{x)  =  6(x)S 
There  are  two  series  which  formally  satisfy  (12),  namely, 

(13)  f(x)  =  -  d(x)  -  0(x+  1)  , 

(13')  /(x)  =  0(x-l)  +  0(x-2)+ 

which,  if  they  converge,  will  be  actual  solutions.  These  will  be  called  the 
direct  sums  to  the  right  and  left  respectively. 

Besides  the  direct  sums,  a  formal  solution  of  (12)  exists  in  the  form  of  a 
contour  integral  due  to  Guichard.f  We  shall  use  it  in  the  form  adopted  by 
Carmichael.J  Let  d  (x)  be  analytic  outside  of  the  region  D ,  extending  to  the 
right  and  lying  within  a  finite  distance  of  the  positive  half  of  the  real  axis 
(see  fig.  1);  then  if  the  contour  L  =  oo  AB  oo  is  drawn  so  as  to  enclose  D, 
and  passes  between  x—  1  and  x  in  the  manner  shown  in  the  figure,  the  function 


//   \       C  Ht)dt 
(14)  f{x)  =  ie^-^-1 


* 


The  solution  written  in  the  form  2  8  ( x )  is  called  a  sum  formula  where  the  operation  2  is 


such  that  %  0  (x)  —  %  0  (x)  =  0  (x)  . 

X+l  X 

f  Wallenberg  und  Guldberg,  loc.  cit.,  p.  12. 
%  Loc.  cit.,  p.  119. 
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formally  satisfies  (12)  by  Cauchy's  formula,  since/(a:+  1)  —  f  (x)  reduces 
to  the  integral  around  ABB'  A' ,  and  the  only  singularity  of  the  integrand 
within  this  region  is  at  t  =  x.* 

If  now  0(0  decreases  along  A  oo  and  B  oo  in  such  a  way  that  the  integral 
converges,  (14)  will  furnish  an  actual  solution  of  (12). 

It  is  obvious  that  in  case  the  region  D  extends  to  the  left,  we  can  take  a 
contour  L  extending  to  the  left  so  as  to  again  enclose  D .  The  sign  of  the 
integral  must  in  this  case  be  changed  if  we  integrate  so  as  to  keep  the  area 
within  L  on  our  left. 


Fig.  1. 

We  therefore  have  four  formal  solutions  of  an  equation  of  the  type  (12), 
the  direct  sums  to  the  right  and  left  respectively,  and  contour  integrals  about 
contours  extending  to  the  right  or  left. 

When  one  of  the  above  sums  is  applied  to  (11)  we  see  that  (9)  will  furnish 
a  formal  solution  of  (1).  Our  main  problem  therefore  is  to  determine  whether  we 
can  so  choose  the  system  of  solutions  faj  (x) ,  Uj  (x) ,  •  ■  •  ,  l„j  (x)  (j  =  1,2, 

•  •  •  ,  n) ,  of  (2),  that  the  application  of  one  or  more  of  the  four  sums  to  (11)  ivill 
lead  to  an  analytic  solution  of  that  equation  and  therefore  of  (1). 

§  3.    The  first  principal  solution  for  /jl  >  0. 

We  choose  first  as  the  solutions  of  (2)  the  first  principal  solutions,  h\j  ( x ) , 
hij  (x) ,  •  ■  •  ,  hnj  (x)  ( j  —  1,2,  •  •  ■  ,  n) .  Take,  further,  as  the  solution  of 
(11)  the  direct  sum  to  the  right.    Thus  we  have  formally 

(10)    u,{x)    (    1)       ^       H(x+1)       +       H(x  +  2)       +  J' 

where  H  (x)  is  the  determinant  |  Ay  (a)  |,  and  Hj  (x)  is  the  minor  of  this 
with  reference  to  the  element  in  the  last  row  and  jih  column. 

*  In  case  6  (  x  )  is  analytic  within  D  ,  the  integral  in  (14)  reduces  to  (13)  by  a  familiar  theorem 
in  the  theory  of  residues,  since  the  only  singularities  of  the  integrand  are  sXx,x-\-l,x-\-2, 

•  •  •  .    That  is,  we  are  led  to  (14)  on  writing  (13)  as  a  series  of  residues. 
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Each  of  the  terras  of  (15)  will  be  analytic  save  at  the  zeros  of  H(x) ,  and 
the  poles  of  b  ( x  +  1 )  and  Hj  ( x  -j- 1 ) ,  or  points  congruent  to  these  points  on 
the  left. 

As  has  been  stated,  hij  ( x )  is  asymptotic  to  s#  ( x )  in  the  right  half  plane, 
so  that  H  (x)  is  asymptotically  represented  by  the  determinant  (x)  |. 
From  (7)  it  follows  that  H  ( x )  will  have  no  zeros  in  the  right  half  plane  at  a 
sufficient  distance  from  the  origin.  Since  the  quantities  ha  (x)  form  a 
system  of  solutions  of  the  homogeneous  system,  we  see  easily  that  H  ( x  +  1 ) 
is  equal  to  A  ( x )  H  ( x ) ,  where  ^4  (a; )  is  the  determinant  |a#  (a;)  |  .  Therefore 
we  have  H  (x  —  1)  =  H  (x) / A  (x  —  1),  and  since  A  (x)  has  no  singu- 
larities in  the  finite  part  of  the  plane  exterior  to  a  circle  of  radius  R ,  it  follows 
that  H  (x)  will  not  vanish  in  the  finite  plane  at  a  sufficient  distance  from  the 
negative  real  axis,  say  at  a  distance  greater  than  R' . 

From  the  situation  of  the  singularities  of  (x)  it  is  apparent  that  Hj  (x) 
is  analytic  in  the  finite  plane  at  a  distance  greater  than  some  quantity  R  from 
the  negative  real  axis,  while  b  ( x  +  1 )  is  analytic  if  x  +  1  lies  outside  a  circle 
of  radius  R  . 

Let  T  be  the  greatest  of  the  quantities  R,  R' ,  and  R.  About  the  origin 
draw  a  circle  of  radius  T ,  and  the  half-lines  tangent  to  it  and  parallel  to  the  neg- 


V 

p 

u 

Fig.  2. 


ative  real  axis  (fig.  2).  Denote  by  P  the  region  enclosed  by  these  lines  and  the 
semi-circumference  in  the  right  half  plane.  Then  if  x  is  any  point  in  the  finite 
plane  exterior  to  P ,  all  the  terms  of  (15)  are  analytic.  If  x  is  within  the  region 
P ,  some  of  the  first  terms  may  have  singularities,  but  all  the  terms  from  the 
rath  remain  analytic,  if  m  be  so  chosen  that  x  +  ra  is  exterior  to  P . 

It  remains  now  to  discuss  the  convergence  of  (15).  Consider  the  rath  term. 
We  will  show  that  it  is  less  in  absolute  value  than  the  mth  term  of  a  rapidly 
converging  series.  To  do  this,  replace  the  elements  of  Hj  ( x  +  ra )  and  H(  x-\-m) 
by  their  asymptotic  forms.  When  this  is  done  it  is  evident  that  we  can  remove 
the  factor  (a;  +  mY(x+m'+rk  ( Pk  )x+m  from  the  kth  column  (k'=  1,  2, 
•  •  •  ,  n)  of  //  (x  +  ra) .    All  these  factors  can  be  removed  from  Hj{x-\-  ra) 
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except  the  one  for  k  =  j ,  so  that  we  have* 

.        b(x+m)Hj(x  +  m)   (x  +  m)v  ( 3  ,  d 

{    }  H{x  +  m)         ~  (x+w)'i(aH^)+^(pie-'i)I+m 

where 

(16')  di  =  h-f, 


\      x+m  ) 


Sj  being  the  minor  of  the  element  in  the  last  row  and  jih  column  of  S  =  |      \ . 

If  we  put  x  +  m  =  mi,  we  see  that  the  typical  term  of  our  series  can  be 
written 

m[Cmi 


mtmi 


M  (mi), 


where  c  and  C  are  constants,  and  M  (mi )  is  a  uniformly  limited  function 
for  m  large  and  x  in  the  vicinity  of  any  point  Xi .  Since  mi  increases  indefi- 
nitely in  the  successive  terms  of  the  series,  it  is  evident  that  the  above 
expression  is  the  general  term  of  a  uniformly  and  rapidly  converging  series, 
irrespective  of  the  values  of  c  and  C ,  provided  only  that  n  >  0 ,  which  we  are 
assuming  for  the  present. 

By  Weierstrass's  theorem  the  series  (15)  will  therefore  represent  an  analytic 
function,  provided  x  is  exterior  to  P;  for  this  makes  all  the  terms  analytic. 
If  x  is  within  P,  the  first  terms  of  (15)  may,  as  seen  above,  have  poles,  the 
remaining  part  of  the  series  converging  uniformly,  and  accordingly  co,  (a;) 
may  have  a  pole. 

The  values  of  coi  ( x ) ,  •  •  •  ,  w«  ( x )  when  used  in  (9)  will  furnish  a  solution 
of  (1)  which  is  analytic  in  the  finite  part  of  the  plane  save  for  poles  within  P . 

The  solution  which  is  obtained  in  this  manner  can  be  extended  to  the  left 
as  follows :  let  Aji  (x)  be  the  minor  of  the  element  in  the  jth  row  and  iih 
column  of  A  (x) ,  and  solve  (1)  for  gi  (x) ,  ■  ■  • ,  gn  (x) ,  as  we  may  do  since 
A{x)  ^  0.    This  changes  the  equations  to  the  form 

(17)   gdx)-±(-iy+i~j^gj(x+l)  +  (-ir^J^b(x+l) 

(i  =  1,  2,  n). 

The  solution  can  now  by  (17)  be  analytically  continued  to  the  left. 

From  this  form  of  (1)  we  can  also  determine  the  location  of  the  singularities 
of  the  solution  without  reference  to  those  of  hij  ( x ) .  Since  the  <7's  are  analytic 
in  the  finite  plane  exterior  to  P ,  it  follows  that,  in  their  continuation  to  the 
left,  singularities  will  be  introduced  only  within  P  at  the  singularities  of 
b  (  x  -f- 1 )  ,  at  the  minors  of  A  ( x  )  ,  and  at  the  zeros  of  A  ( x ) .  A  singularity 
will  in  general  recur  at  points  congruent  on  the  left. 

*  The  series  asymptotically  representing  the  sum,  product,  or  quotient  of  any  number 
of  functions  is  found  by  adding,  multiplying  or  dividing  their  respective  asymptotic  series 
according  to  the  ordinary  rules  for  convergent  series,  Borel,  Les  Series  Divergentes,  p.  30. 
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By  using  the  first  principal  solutions  of  the  homogeneous  system,  and  the  direct 
sum  to  the  right,  we  have  thus  shown  the  existence  of  a  first  principal  solution 
gu  (#)  >  02i  (#) ,  •  •  •  ,  gnl  (x) ,  of  (1),  which  is  analytic  in  the  finite  plane,  save 
for  poles  near  the  negative  real  axis,  congruent  on  the  left  to  the  singularities  of 
an  ( x ) ,  b  ( x  +  1 ),  and  the  zeros  of  A  ( x ) . 

§  4.    The  asymptotic  form  of  the  solution  in  the  right  half  plane. 

We  shall  now  investigate  the  asymptotic  form  of  the  solution  gn  (x) , 
02i  (#)  j  •  •  •  >  gni  (x) .  Manifestly  for  this  purpose  it  is  merely  necessary 
to  consider  more  closely  the  terms  of  (15).  The  asymptotic  form  of  the  typical 
term  has  already  been  determined,  and  is  given  in  (16). 

Before  proceeding  it  will  be  convenient  to  introduce  a  simplification  in 
notation.  When  /  ( x )  is  asymptotically  represented  by  a  power  series  s  ( x ) 
of  negative  integral  powers  of  x,  instead  of  replacing  f  (x)  by  the  expression 
which  is  obtained  on  breaking  off  the  series  s  (x)  at  some  point  and  inserting 
the  remainder  coefficient,  (see  footnote,  §  1),  we  shall  put 

f  (x)  =  [a0;  x], 

giving  in  each  case  the  constant  term  a0  and  the  argument  x  in  the  series.* 


The  leading  factor  of  (16)  can  be  written 


(18) 


x»x  ( pj  tr*  Y  xr'+li       ,     7    ,  my;+*(^> 


the  first  factor  of  which  is  independent  of  m .  We  shall  omit,  for  the  present, 
this  factor,  so  that  we  have  to  consider  the  series 

"  (Rv-  I  p  )m 

(19)  2.  /  T 


m=l  , 

^.n  (to— 1) 


Consider  in  the  first  place  the  first  k  terms  of  (19),  k  <  |  x  | ,  that  is,  the 
finite  series 

i~  '  '  '  "t"  f  i_\  r,—v+n  (z+fc)  \-aJ>  X~\-  K  \  , 

^.M  (fc— 1) 


*  This  notation  is  similar  to  that  used  by  Birkhoff,  these  Transactions,  vol.  9 
(1908),  pp.  373-395. 
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and  let  us  examine  the  first  term  in  detail.    We  have 

2(r,— v)+il 

e*  j  pj  e*-{o-H-*(*+i)}  iog  (1+1/x)     e-     2x  i      xx  Xo 


H)' 


Pj  Pj  Pi  -wx  ■  x< 


=  -  +  —  4-  - 


If  then,  in  the  series    •  +  dj0  /  ( x  +  1 )  +  •  •  •  ,  we  put 

x+l     x\  )■' 
the  first  term  of  (20)  can  be  written 

dj  ,  £1  ,         ,  Lk  +  Erf  ( a; ) 
Pi     x       ' ' 

where  E^k  (x)  is  a  uniformly  bounded  function  for  \x\  sufficiently  large  and  x 
in  the  right  half  plane,  since  it  is  composed  of  remainder  terms  in  uniformly 
convergent  series  and  the  remainder  coefficients  in  asymptotic  series.*  The 
terms  of  (20)  from  the  second  on  can  be  treated  in  the  same  way;  but  we  notice 
that  their  expansions  begin  with  x~*  ,  x~2*  ,  •  •  •  ,  x~<-k~1)lJ'  /  respectively.  The 
first  k  terms  can  therefore  be  written 

il  J_  h  _1_  h  J_  l    hk+Vy-k  (X) 

pj  +  x  +  x2  ~*         h  x*k 

where  rj^k  (x)  is  limited  for  x  large  in  the  right  half  plane. 

It  remains  now  to  consider  the  portion  of  (19)  from  the  term  m  =  k  4-  1 
on.  We  will  find  an  upper  limit  for  the  sum  of  these  terms.  In  the  first 
place  we  notice  that  the  factor 


rj—v+n  (x+m) 


which  occurs  in  the  denominator,  is  greater  in  absolute  value  than  some  quan- 
tity K,  and  in  general  increases  with  m,  since  x  is  in  the  right  half  plane. 
Also  the  expression  [  dj ;  x  +  m  ]  can  be  replaced  by  a  quantity  larger  in  ab- 
solute value.  Consequently  the  portion  of  our  series  from  the  (&+  1  )th 
term  on  is  less  in  absolute  value  than  the  series  whose  general  term  is 


M 


|£  I      |m  (to— 1) 


Pj 


m 

) 


*  To  prove  that  /  (  x  )  ~  a0  +  ai  /  x  +  a2 I  x2  +  •  •  •  ,  we  need  show  only  that  f  *  (  x  )  in 

fix)  =  at,  +  ai  /  x  +  ■  •  ■  +  t  at  +  ft  (  x  )  ]  /  xk    remains    bounded    for   every   k  ,  since 

f  k  (  x  )  =  [ak+i  +  fi+i  ( x )  ]  /  x ,  and  hence  Km  f  4  (  x  )  =  0  .    This  principle  will  be  used  in 

2=00 

what  follows. 
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where  M  does  not  increase  with  \x\  and  K  does  not  decrease.  This  is  a  geo- 
metrical progression,  whose  ratio  can  be  made  as  small  as  we  please  by  taking 
|  a;  |  large.  The  sum  of  the  comparison  series  is  therefore  of  the  order  of  its 
first  term,  which  is  1/  |  x  \*k ,  since  we  are  starting  with  m  equal  to  k  +  1  • 
This  is  of  the  same  order  as  the  last  term  of  the  expression  found  above  for 
the  first  m  terms,  and  consequently  (19)  can  be  written  [dj  /  py,  x] ,  where  the 
series  can  be  taken  to  as  many  terms  as  desired,  since  k  was  arbitrarily  large. 

An  expression  for  Wj{x)  will  be  obtained  by  multiplying  [dj/  pj;  x]  by 
(  —  1  )n+i+1f  (see  (15)),  and  by  the  first  factor  in  (18),  which  was  removed  from 
(19). 

The  quantity  we  wish  to  determine  is  coj  (x)  ha  ( x ) .    Therefore,  since 

hij  (x)  =  x*x(pj  e~'1  Y  xr>  [  Sij;  x  ] ,  we  have 

(21)  Wj  ( x )  hij  ( x )  =  aT-*  [  (  -  1  )n+j+1  dj  Sij I  Pj;x]. 
Letting  j  =  1,2,  •  •  •  ,  n  and  adding  one  finds 

(22)  (7a  (x)  =  xv-»[hi;  x]  «  =  1,  2,  •  ■  - ,  n) . 
where 

\     Pi  P2  Pn  ) 

From  (6)  and  (16')  we  obtain 

(23)  h  =  — —  X)  (  —  1       Pi  P2  •  •  •  Pj-i  Pj+i  ■  ■•  Pn  Sij  Sj 

( i  —  1 ,  2  ,  •  •  •  ,  n  )  . 

Let  us  now  try  to  find  series  which  will  formally  satisfy  (1).  Write 


+  '••)       «  =  1>  2-  n)' 


(24)  ^(ix)  =  a: 

and  substitute  in  (1).  When  we  develop  the  resulting  equations  according  to 
decreasing  powers  of  x  it  is  evident,  on  remembering  that  p.  >  0 ,  that  we  can 
uniquely  determine  the  successive  constants  gif  ,  ■  •  ■  by  comparing  coef- 
ficients. The  developments  in  all  the  equations  start  with  x";  equating  its 
coefficients  we  find 

n 

0  =  E«i3  91  ( z  =  1 ,  2  ,  •  •  •,  n  -  1 ) , 

i=i 

From  this  we  have  for  the  first  coefficient  in  <jrt-  ( x ) 

hA 

(25)  0i=  (~  l)w+m^  «-l,2,  '",»). 
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where  Am  is  the  minor  of  the  determinant  A  =  |  ai;- 1  with  reference  to  the  ele- 
ment in  the  last  row  and  tth  column. 

We  call  the  series  gi  (x) ,  •  •  •  ,  gn  (x)  the  formal  solution  of  (1).  These 
series  are  obtained  by  direct  substitution,  and  in  general  diverge. 

Since  the  expansions  in  the  expression  for  gn  (x)  and  in  the  formal  series 
g%{x)  ( i  =  1,2,  •  ■  •  ,  n )  are  of  the  same  form,  the  question  naturally  arises 
as  to  whether  their  coefficients  are  the  same.  We  can  show  that  this  is  the 
case  in  the  following  way.  We  know  that  gn(x) ,  •  •  •  ,  gnl  ( x )  as  given  in  (22) 
satisfy  (1),  and  if  we  substitute  this  form  of  the  solutions  in  (1)  and  compare 
coefficients  of  quantities  of  the  same  order  we  see  that  we  must  have  precisely 
the  same  equations  satisfied  by  h{ ,  h^,  ■  ■  ■  that  were  satisfied  by  g{ ,  g[l) ,  •  •  • 
(i  =  1,2,  •  •  •  ,  n) .  Since  these  equations  had  unique  solutions,  this  means 
that  corresponding  coefficients  in  the  series  gn  ( x )  and  gi  (x)  (i  —  1 ,  2 ,  •  •  •  ,  n) 
are  equal. 

Thus  we  see  that  the  first  principal  solution  is  asymptotically  represented 
by  the  formal  solution  in  the  right  half  plane. 

Although  we  have  shown  easily  that  the  constants  in  the  series  for  gn  (x) 
must  be  the  same  as  those  in  the  formal  series  gi  ( x ) ,  it  is  interesting  to  note 
that  we  can  show  by  direct  computation  that  hi  is  identical  with  gi .  To  do 
this  we  must  eliminate  the  quantities  pj  and  Sij  from  (23). 

When  we  substitute  the  formal  series  (3)  in  (2)  and  compare  coefficients  of 

xU.+r]  we  £n(J 

n 

X!  ( Oit  -  Sik  pj)  skj  =  0  (i  =  l,2,  •■•  ,  n), 

k=l 

where  So,  =  0  if  i  4=  k ,  and  8  a  =  1 .  Since  pj  is  a  solution  of  (4),  the  last  equa- 
tions have  solutions  other  than  8hj  —  0 .    The  first  n  —  1  of  them  give 

Pi  Sij  =  ^laikSkj  ( *  -  1  >  2,  •  •  ■  ,  n  -  1 )]. 

If  now  we  introduce  the  quantities  pi ,  •  •  •  ,  pj-i ,  pj+i ,  •  •  •  ,  pn  into  the  succes- 
sive columns  of  the  determinant  Snj  and  eliminate  the  p's  by  means  of  the 
above  expressions,  we  obtain  the  relation 

(26)       Pl  •  •  •  Pj-i  Pj+l  •  '  ■  Pn  Snj  =  A,a  Sij  +  An2  S23  +  ■  "  '  +  Ann  Snj, 

where  Ay  and  Sij  are  the  minors  of  A  and  S  respectively  with  reference  to  the 
element  in  the  ith  row  and  jth  column,  and  where  we  have  put  Snj  =  Sj, 
the  two  being  identical  by  definition. 

If  in  (26)  we  let  j  =  1,2,  ■  ■  ■  ,  n,  and  substitute  in  (23),  we  find  for  i  =  1 

h  =  — ^2 —  2^     ^  ( —  1 )    *w     =  ~2  =  gi ' 
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since 

£T      '     1J*kj     10  if  k  4=  1. 

In  a  similar  way  we  can  show  hi  =  <7»  (i  =  2 ,  3,  •  •  • ,  n) . 

We  summarize  the  results  of  this  and  the  preceding  section  in  the  following : 

Theorem  I.  If  fx  >  0 ,  the  system  of  non-homogeneous  linear  difference  equa- 
tions (1)  has  a  first  principal  solution  gn  (x) ,  g2i  (x) ,  •  •  •  ,  gn\  (x) .  All  the 
constituent  elements  of  this  solution  are  analytic  throughout  the  finite  plane,  save  for 
poles  at  the  poles  of  the  functions  aij  (x)  ,b  (x  +  1 ) ,  at  the  zeros  of  the  determi- 
nant |  a^  ( x )  | ,  and  at  points  congruent  to  these  points  on  the  left.  There  exist 
n  series  gi  (x) ,  g2  (x) ,  ■  ■  •  ,  gn  (x)  which  formally  satisfy  the  system  (1),  and 
whose  coefficients  can  be  determined  by  direct  substitution.  These  series  in  gen- 
eral diverge,  but  when  x  approaches  infinity  in  any  right  half  plane,  ga  (x)  is 
represented  asymptotically  by  gi  ( x ) . 

This  theorem  is  supplemented  by  Theorem  II  of  §  6. 

§  5.    A  formal  property  of  the  direct  sum. 

The  solution  of  (1)  which  we  have  found  was  determined  by  making  use  of 
one  of  the  sets  of  principal  solutions  of  (2).  The  question  now  arises  as  to 
whether  a  different  solution  of  (1)  would  be  obtained  if  we  used  a  different 
fundamental  system  of  solutions  of  (2)  and  employed  again  the  direct  sum 
to  the  right.    The  answer  is  in  the  negative  as  will  be  shown. 

Let  hj  ( x ) ,  l2j  ( x ) ,  •  •  •  ,  lnj  (x)  (j  =  1,  2,  •  •  •  ,  n)  be  any  other  funda- 
mental system  of  solutions  of  (2).  From  the  theory  of  homogeneous  systems 
we  have 

ii 

hi  {%)  =  Z!  <Pij  ( x )  hu  ( x ) , 

i=l 
n 

(27)  hj(x)  =  J2<Pij(z)h2i(x),  U  - 1,  2,  n), 

i=\ 


lnj(x)  =  X  <Pij  (  X  )  hni  (X), 
i=l 

where  the  functions  <pij  ( x )  are  periodic  of  period  1 ,  and  the  determinant 
$  (x)  =  |  <pi3-(x)  |  +  0. 

Denote  by  gi  (x)  (i  =  1,  2,  •  •  • ,  n)  the  solution  of  (1)  arising  from  (27)  . 
Then  one  obtains  the  formula  analogous  to  (9), 

n 

(28)  gi(x)  =  T,6j(x)lij(x)  (i  =  l,2,  ...,rO, 
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where  the  functions  03-  (x)  ( j  —  1 ,  2 ,  •  ■  •  ,  n)  are  to  be  determined.  As 
before,  let  L  (x)  denote  the  determinant  |  Uj  (x)  | ,  and  Lj  ( x )  the  minor  of 
this  with  reference  to  the  element  in  the  last  row  and  jth  column;  then  we 
have  by  (15),  since  the  work  to  that  point  was  purely  formal, 

(29)    °A*)=(-1Y^[HX+L\1L;(1]+11+---]    0-1,2,  -.-). 

From  (27)  it  follows  that  L  ( x )  =  H  ( x )  •  $  ( x ) .  Similarly,  we  can  show 
that 

Lj  (x)  =  *iy  (x)  H,  (x)  +  $2i  (x)  H2  (x)  +  •  •  •  +       (a.')  ff„  (z) 

(;  =  1  ,  2,  ■■■  ,  n), 

where  ( a; )  is  the  minor  of  <3?  ( x )  with  regard  to  the  element  in  the  ith  row 
and  jih  column,  and  Hj  (x)  is  the  function  in  §  3.  Therefore,  since  |  <3?  (x)  | 
and  all  its  minors  are  periodic  functions,  we  obtain 

(a)  =     *(*).    §*«  {X){  ,H{z+l)  +  •  •  • )    (J  =  1>2'---'n)' 

or  by  virtue  of  (15), 

Substitute  the  values  just  obtained  for  6\  (x) ,  •  •  •  ,  6n  (x)  in  (28);  we  find 

+  (-l)-1^^g(-l)*fl**(a)i«(*>. 

If  we  make  use  now  of  the  values  of  (x)  as  given  by  (27),  the  last  expression 
becomes,  upon  rearranging, 

?f(*)-fl?rt*«(a!)t(-l)ww»(*)«i*(ai) 

\il  )   j=l  fc=l 
+  


+  <-  U^t^E^CarjEC-  l)^^(-r)*nA-(.T). 


224  K.  P.  WILLIAMS:  [April 

Since 

±(-  »*»(.)•-(.)  =  {(1^™*^),  if  jm ,, 

the  above  expression  for  </;  ( x )  easily  simplifies,  and  finally  we  have 

</;  (a)  =  COl  (a)  /in  O)  +  C02  (JC)  A*2  O)  +  •  •  •  +  Wn  (x)  hin  (x)  . 

Therefore  ~gi  ( x )  =  ^1  ( a; )  ( i  =  1 ,  2 ,  •  •  •  ,  n ) ,  so  that  we  have  the  same 
solution  obtained  in  §  3. 

This  result  can  be  stated  as  follows:  When  the  equation  (11)  can  be  solved 
by  means  of  a  direct  sum  to  the  right,  the  solution  of  (1)  obtained  by  the  method  of 
variation  of  constants  is  independent  of  the  choice  of  the  fundamental  system  of 
solutions  of  the  associated  homogeneous  system,  for  in  this  case  the  series  which 
occur  are  convergent,  and  therefore  can  be  rearranged  as  above. 

§  6.    The  asymptotic  form  of  the  solution  in  the  left  half  plane. 

By  an  application  of  the  results  of  the  last  section  we  can  determine  the 
asymptotic  form  of  the  first  principal  solution  in  the  left  half  plane. 

Take  as  the  solutions  of  (2)  the  intermediate  solutions  h\j  ( x ) ,  h'2j  ( x ) ,  •  •  • , 
Ki(x)  (j=  1,2,  •••,«),  of  §1.  Put 

n 

(30)  g'n(x)  =  Y,aj(x)h'iJ(x)  (t-1,2,  •..,»), 

3=1 

and  apply  the  direct  sum  to  the  right  in  order  to  determine  co[  ( x ) ,  •  •  •  , 
Z)'n  (x);  a  formal  value  of      ( x )  will  be 

(si)       ;;WM-W6<'t1^(*+1>  +  ...1 

where  IV  (x)  =  |  h'ij  (x)  \,  and  H'j  (x)  is  the  minor  of  H'  (x)  with  respect 
to  the  element  in  the  last  row  and  jth  column.  The  same  convergence  proof 
used  in  §  3  applies  here,  provided  x  is  sufficiently  distant  from  the  real  axis, 
since  we  have  Ay  ( x )  ~  Sij  ( x )  with  regard  to  v  in  the  right  half  plane. 

According  to  the  result  of  the  last  section  this  new  solution  must  be  identical 
with  the  first  principal  solution.  As  derived  above  it  is  known  to  exist  only 
at  a  sufficient  distance  from  the  whole  extent  of  the  real  axis. 

The  typical  term  of  (31)  will  again  have  the  asymptotic  form  (16),  the  re- 
presentation holding  with  regard  to  x  in  the  left,  and  with  regard  to  v  in  the 
right  half  plane.  To  find  the  asymptotic  form  of  the  series,  divide  it  into  two 
parts :  the  first  k  terms,  where  k  is  such  that  x  4-  k  +  1  =  Xi  is  a  point  in  the 
right  half  plane  near  the  imaginary  axis,  and  secondly,  the  remaining  terms. 

Take  first  the  case  where  x  approaches  infinity  along  any  ray  making  an 
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acute  angle  with  the  negative  real  axis.    The  first  k  terms  can  again  be  written 


in  the  form  (20). 
that  the  quantity 


By  the  definition  of  k  we  have  in  this  case  k  <  |  x  | ,  so 


>',—  H-m  (x+m) 


( m  <  k ) 


which  appears  in  the  denominator  of  the  typical  term  can  be  developed  in  the 
same  way  that  the  first  term  of  (20)  was  treated,  since  all  the  expansions  are 
valid  if  m  <  |  x  \ .  Thus  we  have  immediately  that  the  first  k  terms  of  (31) 
when  multiplied  by  h\j  (x)  can  be  put  in  the  form  (21). 

Consider  now  the  terms  of  (31)  from  the  (k  +  1  )th  on.  They  can  be 
written 

^  b  ( Xi  +  m )  H'j  ( Xi  +  m ) 
m=0  H'  (xi  +  m ) 

Since  Xi ,  X\  +  1 ,  •  •  •  are  points  in  the  right  half  plane,  the  asymptotic  form 
of  the  typical  term  of  this  series  is  again  given  by  (16),  the  representation 
being  now  with  regard  to  v .    The  series  therefore  takes  the  form 


(32) 


(e"/pyriH 


■  Mm  (zi+  m), 


where  Mm  ( a*i  +  m )  remains  uniformly  bounded  for  |  v  |  large. 

Let  x\  =  Mi  +  v  V  —  1,  rj  =  £  +  f  V  —  1;  then  we  can  so  choose  X\  that 
MMi  +  £  —  v  >  0 .  We  can  also  choose  x  so  that  the  sign  of  fiv  +  f  is  the 
same  as  that  of  a .  It  can  then  easily  be  shown  that  the  absolute  value  of  the 
denominator  of  the  (m  +.  1  )th  term  of  (32)  is  greater  than  \Xl^Zl+m')+r~v  \ . 
The  absolute  value  of  (32)  is  therefore  less  than 

m\  (^/p.-riv 


Xl 


M  *!  +  >', 


•It  Pj 


where  M  is  a  constant  not  depending  on  x .  Since  |  e"  /  x  f  p;-  |  can  be  made  less 
than  unity,  the  sum  of  the  last  series  is 


(33) 


\Xi 


where  M'  is  bounded. 

The  terms  of  (31)  from  the  (&+  1  )th  on  are  also  to  be  multiplied 
by   h'ij  ( x ) ,    the    absolute   value   of    the    principal   part   of    which  is 


x 


(Pie-*) 


~Ui\x  |_— (»u  +  f)  argx 


From  (33)  the  absolute  value  of  the  product 


is  seen  to  be  less  than  M'  \  x  |M*+£  |  pje 


—H   Hi— Mi 


This  vanishes  more  rapidly 


than  any  power  of  1  /  x  along  the  rays  considered,  inasmuch  as  |  p} ■  e  u  |  is 
constant;  for  along  these  rays  u  increases  negatively.    Consequently  the  terms 
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of  (31)  from  the  ( k  +  1  )th  do  not  affect  the  asymptotic  form  of  c3J-  (z)  h'ij(x) , 
which  is  therefore  given  by  (21).* 

It  remains  now  to  consider  the  case  where  x  approaches  infinity  parallel  to  the 
negative  real  axis  at  a  sufficient  distance  from  it.  A  slight  modification  in  the 
treatment  of  the  first  k  terms  must  be  made,  for  the  quantity  ( l-\-m/ xYJ~"+>i(x+m) 
cannot,  in  all  the  terms,  be  developed  as  before,  for  in  the  last  terms  \mjx\ 
will  approach  very  close  to  unity.  The  modification  to  be  made  is  apparent: 
we  divide  the  first  k  terms  into  two  parts;  the  first  I  of  them,  where,  for  instance, 
I  is  the  smallest  integer  such  that  |  x  +  I  |  <  §  |  x  |  (such  a  choice  is  always 
possible  when  x  goes  to  infinity  parallel  to  the  negative  real  axis ),  and  secondly 
the  remaining  terms. 

The  first  I  terms  when  multiplied  by  h'ij  (x)  will  again  have  the  asymptotic 
form  (21).  We  can  show  readily  that  the  second  group  of  terms  when  mul- 
tiplied by  the  same  quantity  will  contribute  nothing  to  the  asymptotic  from 
of  co;  ( x )  h'ij  ( x ) .  The  details  of  the  argument  are  similar  to  the  treatment 
of  the  terms  after  the  kth,  and  we  will  merely  outline  them  in  a  note.f 

For  this  case  we  therefore  have  again  that  coy  (x)  h'ij  (x)  is  given  by  (21). 
From  this  it  follows  immediately  that  g'n  (x)  ~  g{  (x)  in  the  left  half  plane. 
Remembering  that  the  new  solution  is  identical  with  the  first  principal  solution, 
we  have: 

Theorem  II.  If  /x  >  0 ,  the  first  principal  solution  is  asymptotic  to  the  formal 
series  in  the  left  half  plane  at  a  sufficient  distance  from  the  real  axis.% 

*  If  f(x)=<p(x)  +  d(x),  where  <p  (  x  )  ~  xn  (  a0  +  ajx  +  •  •  • )  along  some  ray, 
while  6  (x  )  becomes  small  of  order  1  /  a\x  I  (  a  >  1  )  ,  then  /  (  x  )  ~  xn  (  a0  +  ax  /  x  +  •  •  •) 
since  1  /  a1  :t  •  is  smaller  than  any  power  of  1  /  x  for  x  large. 

t  The  terms  in  question  can  be  written  in  the  form  (32),  I  =  m  =  k ,  and  x  replacing 
xi .  We.  can  show  first  that  the  absolute  value  of  the  typical  term  would  be  increased  if  we 
changed  (x  +  m)ixX+rJ-1'  in  the  denominator  to  (x  +  I  )f *+'■/-",  and  the  factor  (x  -\-m)»-m 
to  x*m,  where  xo  is  the  smallest  in  absolute  value  of  x  + 1,  ••• ,  x+k  .  We  obtain  then  for  a 
series  greater  term  by  term  in  absolute  value,  a  geometrical  progression,  which  can  be  summed, 
the  sum  being  less  than 

M\(e*l  p})*+l\ 


since 


|  (x  +  l)^x^rj-"\  ■  \xo\i*1 


<  M\(xl2)"-»-x-rj(e*  I  Pi)x\ 


<  1  ,  and      +     <  |  a;  /  2  | .    The  product  of  this  and  the  absolute  value  of  the 


leading  factor  of  the  asymptotic  form  of  hjj  (a;)  is  easily  seen  to  be 

M\x\v2^u+t 


M\{xl2)v2*x+ri\  = 


2" 


(x  =  u  +  v  V —  1)  .  Since  u  is  practically  equal  to  —  \x\,  and  M  is  a  fixed  number,  this 
decreases  more  rapidly  than  any  power  of  1  /  x  . 

t  It  should  be  noted  that  the  asymptotic  character  of  this  solution  of  the  non-homogeneous 
equation  is  somewhat  different  from  that  of  the  solutions  of  a  homogeneous  equation.  For 
instance,  the  asymptotic  form  of  r  ( x  )  holds  only  for  —  7r  <  arg  x  <  w  .    That  functions 
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It  is  interesting  here  to  refer  to  (17).  The  quantities  A}i  (x)  /  A  (x) 
(j,  %  =  1,  2,  •  •  •  ,  n) ,  which  occur  there  have  zeros  of  order  \x  at  infinity. 
Consequently,  if  we  assume  that  there  exists  a  constant  X  such  that  lim  xK  gi(  x ) 
=  hi,  ki  finite,  as  x  approaches  infinity  parallel  to  the  negative  real  axis, 
equation  (17)  gives 

An  ( ,  ,  ciij  \  x'gjjx) 

hm  — 


Hm  x%  ( x )  =  t  (  -  1  )>+1  ^  (  1  +  |  +  •  •  •  ) 


+  lim  (- 
t  we  must  have 
X  =  fi  -  v ,       hi  =  (-  1 )"+;+! 


From  this  relation  we  see  that  we  must  have 

bAni 


so  that 

9il  (X)  =  Xv~»[(-  l)n+l+lhAf+  M*)]    (i  =  1,  2,  n), 

where  (a')  tends  to  zero  as  |  .r  |  becomes  infinite.  From  (24)  and  (25)  it 
is  clear  this  is  obtained  by  breaking  off  the  series     ( x )  after  the  first  term. 

§  7.    The  second  ■principal  solution  for  fx  >  0 . 

We  shall  now  show  the  existence  of  a  solution  of  (1)  defined  throughout 
the  finite  plane,  and  whose  asymptotic  form  is  known  in  the  whole  left  half 
plane.  It  will  of  course  differ  from  the  first  principal  solution  merely  by  a 
solution  of  the  homogeneous  system  (2). 

Take  as  solutions  of  (2)  the  principal  solutions  on  the  left,  that  is,  hy  ( x ) , 
( x )  j  •  •  •  ,  hnj  ( x )  (j  =  1 ,  2 ,  •  ■  • ,  n )  ,  where  (x)  ~  ( x )  in  the  left 
half  plane.  Put 

n 

(34)  gi2  (x)  =  T,aj(x)hij  (x)  (i  =  l,  2,  n). 

We  find  in  precisely  the  same  way  as  before 

.b(x+l)flj(x+l) 


(35)  Aoj(x)  -  (-  l)n+j 


H  (x+1) 


where  H  (x)  is  the  determinant  |  hij  (x)  \,  while  Hj  (x)  is  the  minor  of  this 
with  respect  to  the  element  in  the  last  row  and  jth  column. 

If  we  attempt  to  solve  (35)  by  the  direct  sum  to  the  left,  which  we  would 
naturally  try  in  this  case,  the  series  obtained  is  found  to  diverge.    We  can, 

actually  exist  with  the  property  specified  in  Theorem  II  can  be  shown  by  considering  the 
series  l/a;  +  l/  (a;+l)2  +  l/(-c  +  2)3+---.  It  is  easy  to  prove  that  this  function  is 
asymptotic  to  its  formal  expansion  in  powers  of  1  /  x  ,  with  regard  to  x ,  for  —  tt  <  arg  x  <  it  , 
and  also  for  arg  x  =  ±  ir ,  provided  that  \v  |  >  a  >  1 ,  where  x  =  u 
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however,  apply  the  integral  formula  (14).  To  prove  this,  we  note  first  that 
the  situation  of  the  singularities  of  Hj  ( x ) ,  and  the  fact  that  H  ( x )  ~  |  (  x )  | 
with  regard  to  v  in  the  right  half  plane,  show  that  the  expression  on  the  right 
in  (35)  is  analytic  outside  some  region  D,  such  as  was  described  before. 
Accordingly  the  expression 


(36) 


5.(jB)aB  (-!).«  f»(«+l)g/(*+l)  *  


will  formally  satisfy  (35),  and  in  case  the  integral  converges  suitably  it  will 
represent  an  analytic  solution.  The  contour  L  is  the  contour  <x>  AB  <x>  used 
earlier. 

To  examine  the  nature  of  the  integrand  we  make  use  of  the  intermediate 
solutions;  for  although,  as  noted  above,  H  (t)  is  asymptotic  to  the  deter- 
minant |  Sij  (t)  \  with  regard  to  v  in  the  right  half  plane,  the  minor  Hj  (t) 
is  not  asymptotic  to  the  corresponding  minor  of  |  s#  ( t)  \ ,  unless  j  =  n, 
(see  §  1).  We  therefore  replace  the  elements  of  Hj  (x)  by  their  values  as 
given  by  (8).  When  we  do  this,  the  resulting  determinant  can  be  simplified  so 
as  to  give  finally 

(37)    Hj  (i)  =  Hj  (x)  +  9m  (x)  H'J+1  (x)  +•••+*>«  (x)  H'n  (x) , 

where,  as  in  the  preceding  section,  H]  ( x )  is  the  minor  of  H'  ( x )  =  |  ( x )  | 
with  respect  to  the  element  in  the  last  row  and  jth  column,  while  <pj+i  (x) , 
•  •  •  ,  (fn  (x)  are  periodic  functions,  and  are  analytic  at  a  sufficient  distance 
from  the  real  axis. 

The  first  factor  of  the  integrand  in  (36)  can  now  be  written 

bjt+DHjjt+l)     A  fr(*+l)gj(*+P 
(dh)  - —  "  =  2-f^i(U —  - —  -,       <pj  (t)  —  I. 

H(t+1)  7^3  H(t+1) 

When  we  make  use  of  the  asymptotic  form  of  the  quantities  involved,  the 
first  term  of  (38)  becomes,  (see  (16)), 

/on\  (e»/Pi)m      (  ,    |  W  ,  Vj(t)\ 


where  r)j  (t)  is  limited,  decreasing  for  large  \  v\  (t  =  u  +  v  V —  1).  When 
|  v  |  is  sufficiently  great  it  is  evident  that  (39)  will  decrease  along  A°o  and  5oo 
in  approximately  the  manner  of  an  exponential  function,  inasmuch  as  we  are 
assuming  p,  >  0 .  Since  the  functions  <pj+i  ( t ) ,  ■  •  ■  ,  <pn  ( t)  are  periodic  and 
limited  along  any  line  parallel  to  the  real  axis,  like  results  hold  for  all  the  terms 
of  (38).  The  quantity  (e2n(x~t)J~1  —  1  )_1  remains  finite  along  L;  conse- 
quently the  integral  will  converge  in  such  a  manner  that  co,  ( x )  is  an  analytic 
function  outside  of  the  region  D . 
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With  these  values  of  (Oj  (x) ,  equation  (34)  gives  a  solution  of  (1)  analytic 
on  the  left  outside  of  D  .  This  solution  can  be  extended  to  the  right  by  means 
of  the  equations  (1)  themselves.  From  the  fact  that  a#  (.r)  and  b  {x  +  1) 
are  analytic  in  the  finite  plane  outside  a  certain  circle,  the  extended  solution 
can  have  singularities  only  at  points  congruent  on  the  right  to  the  singularities 
of  ciij  ( x )  and  b  (  z  +  1 )  • 

We  now  come  to  the  asymptotic  form  of  the  solution.  For  this  purpose  we 
break  up  the  contour  L  into  two  contours,  Li  and  Z2  in  the  way  that  Birkhoff 
has  done.*  Here  L\  is  a  fixed  contour  ocAB<x>,  outside  of  which  the  integrand 
in  (36)  is  analytic,  and  in  case  x  lies  between  A°o  and  £oo  extended  to  the 
left,  Li  is  a  loop  circuit  around  x  and  the  first  point,  x  4-  /,  to  the  left  of  AB , 
that  is  congruent  to  x;  while  if  x  is  above  A°o  or  below  Boo  ,  Z2  is  a  loop  circuit 
to  infinity  parallel  to  the  positive  real  axis.  In  either  case  L\  and  Z2  are 
together  equivalent  to  L ,  and  L\  can  be  deformed  so  that  no  point  congruent 
to  x  falls  on  it. 


■  x 


A 


X  T  x  +  l 


B  00 
Fig.  3. 

Denote  the  parts  of  cij  (x)  arising  from  Li  and  X2  by  (x)  and  co;")  (x) 
respectively.  Then 

(40)  gi2  ( x )  =  Z  T  ( x )  ha  ( .r )  +  Z  %P  (x)hij(x)    «  -  1 ,  2 ,  •  ■  • ,  n) . 
Consider  first  the  case  where  x  lies  between  A<x>  and  .Boo;  we  have 

(41)  -fix)  =  (-  iy^$j»(«+»>g>(«+->. 

m=i        jH  ( a;  -+-  m ; 

The  asymptotic  form  of  the  typical  term  is  again  given  by  (16).    Also  we  see 
*  The  reason  for  this  is  evident.    The  singularities  of  the  integrand  in  (36)  are  at  the  con- 
gruent points  x,x  +  l,x  +  2,---,  due  to  the  zeros  of  e^(x-t)S-i  _  ^  ancj  ^e  Qxq^  singu- 
larities of  k  ^  ^       )      ( ^  "l~  ^  )   which  are  in  the  vicinity  of  the  positive  real  axis.    The  resi- 

dues  at  the  first  points  will  produce  a  series  of  the  form  (13),  and  the  sum  of  the  residues  at 
the  last  points  will  be  a  periodic  function.  It  is  therefore  natural  to  consider  the  two  different 
groups  of  terms  separately. 
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that  precisely  the  treatment  applied  to  the  first  k  terms  of  (31)  for  the  case 
where  x  approaches  infinity  parallel  to  the  negative  real  axis,  holds  for  (41), 
provided  AB  is  taken  far  enough  to  the  left.  Therefore  Zip  ( x )  hij  (x)  is 
again  given  by  (21),  so  that  we  have 


n 


(42)  £  wf(x)hij  (z-)~gi(x)  (*  =  l,  2,  •■■ ,  n), 


3 

3=1 


for  x  between  Ax>  and  Bx>. 

Next  let  x  lie  above  Aoo  or  below  Boo;  then  (41)  is  replaced  by  the  infinite 
series 

gp{x)  =  ( _  x  }»+m  ±  6(*  +  m)gj(*+m) 

™=i  i7(  .r  +  m ) 

Since  the  fundamental  system  hij  (x) ,  h2j  (x) ,  ■  ■  ■  ,  hnj  ( x )  can  be  expressed 
in  terms  of  the  solutions  hij  (x) ,  h2j  (x) ,  •  •  ■  ,  hnj  (x)  (j  =  1,2,  •  •  •  ,  n) , 
it  follows  from  the  results  of  §  5  that 

n  n 

(43)  Zwf  (a;)M»0  =  £<»i(a:)MaO  =  9a  (x)     (t"  =  x>  2'  n>- 

But  we  have  found  that  gn  (x)  ~  gf»  (.r)  along  any  ray  making  an  angle 
with  the  negative  real  axis,  or  parallel  to  that  axis  at  a  sufficient  distance 
from  it.    We  see  then  that  (42)  is  true  throughout  the  left  half  plane. 

We  still  have  to  examine  the  quantity  ufp  (x) .  L\  is  a  fixed  contour,  so 
that  co^1'  (a;)  is  periodic,  and,  moreover,  its  absolute  value  will  have  a  finite 
upper  limit.    Furthermore,  we  can  write 

(-  i)^  r  b(t+i)Hj(t+i)  _dt 

The  integral  approaches  a  constant,  Cj ,  when  x  approaches  infinity  along  the 
negative  imaginary  axis,  so  that 

(44)  lim  e2^"1  fy»  (  X  )  =  Cj  (  for  arg  x  =  ~  )  • 
From  (36)  itself  we  see 

(45)  lim  up  (x)  =  Cj  ( for  arg  x  =  ^  J  , 

.r=oo 

where  Cj  is  constant.* 

Consider  now  the  product  ufp  (x)  (x) .  It  can  easily  be  shown  that 
the  leading  factor  x*x  {pje~'1  )x xr>  in  the  asymptotic  form  of  hij  (x)  vanishes 
more  rapidly  than  any  power  of  1  /  x  along  any  ray  in  the  left  half  plane,  includ- 
ing the  positive  axis  of  imaginaries.  Then,  since  co^0  (x)  remains  bounded 
(x)  hij  (x)  will  also  vanish  more  rapidly  than  any  power  of  1  /  x ,  and  con- 

*  The  inversion  of  order  of  the  double  limits  employed  here  is  easily  justified. 
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sequently  it  will  not  effect  the  asymptotic  form  of  a>;  ( x )      ( x ) ,  so  that 

gi%(x)  ~  gi{x)  (!^argz<y). 

By  using  the  second  principal  solutions  of  the  associated  homogeneous 
system,  and  expressing  the  sum  formulae  by  contour  integrals  to  the  right,  we 
have  thus  proved: 

Theorem  III.  If  /x  >  0 ,  there  exists  for  (1)  a  second  principal  solution  gn  ( x ) , 
#22  (x)  >  •  •  •  >  Qni  (x)  analytic  throughout  the  finite  plane,  save  for  poles  con- 
gruent on  the  right  to  the  poles  of  an  (x)  and  b  (x  +  1).  This  solution  is 
such  that  #i2  (x)  ~  gi  (x)  when  x  approaches  infinity  along  any  ray  in  the  left 
half  plane,  inclusive  of  the  positive  imaginary  axis. 


§  8.    The  relation  between  the  principal  solutions. 

The  difference  between  any  two  solutions  of  (1)  is  a  solution  of  the  homoge- 
neous system  (2).  To  determine  the  relation  between  the  principal  solutions 
we  make  use  again  of  the  second  principal  solutions  of  the  homogeneous 
system.    We  know  that  we  must  have 

n 

ga  (x)  —  ga  (x)  =  53  Qj  (x)  ha  (x)      (2  =  1,2,  ■••,»), 

where  the  functions  6i  ( x ) ,  •  •  •  ,  6n  ( x )  are  periodic.  We  have  really  already 
determined  these  periodic  functions  in  the  form  of  contour  integrals;  for  a 
comparison  of  (40)  and  (43)  with  the  above  expression  shows  immediately  that 
we  must  have  6j  (x)  —  co^11  (x)  when  x  is  sufficiently  distant  from  the  real  axis. 

The  functions  6i  ( x ) ,  •  •  •  ,  6n  ( x )  are  therefore  analytic  at  a  sufficient 
distance  from  the  real  axis,  and  as  they  have  no  essential  singularities  in  a 
period  strip,  they  are  rational  functions  of  eUx  Thus  we  can  write 

mj    Ik  7)00   

where  wij ,  rij,  nj ,  and  are  positive  integers,  and  Dm  ,  are  constants, 
while  akj  is  a  fixed  point. 

When  we  remember  the  identity  of  dj  ( x )  and  Z>1} }  ( x ) ,  for  x  sufficiently 
distant  from  the  real  axis,  we  see  from  (44)  that  we  must  have 

df  =  0 ,       h  ^  0 

,.A,  (.7  =  1,2,  n). 

(46) 


*  Osgood:  Lehrbuch  der  Funktionentheorie,  I,  1st  edition,  p.  404,  or  2d  edition,  p.  466. 
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df=0,  k<0, 


(47) 


U  =  1,  2, 


n). 


=  Cj 


The  relation  between  the  principal  solutions  is  accordingly  given  by 


(48) 


(Jiz  ( x )  -  ga  (x)  =  XI  hj  ( £ )  £  £ 


n    _  m,-  qk 


j=i  k=i  i=i  ie 


i  =  1 ,  2  ,  •  •  •  , 


rc), 


with  the  conditions  (46)  and  (47). 


9.    The  case  m  <  0 . 


We  have  thus  far  explicitly  assumed  that  p.  is  positive.  We  shall  now  in- 
dicate the  treatment  for  ix  negative. 

The  series  (15)  diverges  in  this  case,  and  so  the  direct  sum  to  the  right  no 
longer  furnishes  a  solution;  but  if  we  apply  the  direct  sum  to  the  left  to 
(35),  we  can  show,  by  using  the  asymptotic  form  of  the  terms,  that  the  series 
so  obtained  will  converge  uniformly  at  a  sufficient  distance  from  the  positive 
real  axis. 

By  means  of  a  direct  sum  to  the  left  we  thus  establish  the  existence  of  a 
solution  analytic  in  the  finite  plane  save  for  poles  at  the  poles  of  a#  (x)  and 
b  ( x  +  1 )  and  points  congruent  on  the  right.  This  solution  corresponds  to 
what  was  called  the  second  principal  solution  when  n  is  positive. 

If  now  we  take  the  first  principal  solutions  as  solutions  of  the  homogeneous 
system,  it  can  easily  be  shown  that  the  resulting  sum  formula?  can  be  taken 
to  be  contour  integrals  about  a  contour  to  the  left.  This  gives  a  second 
solution  of  the  non-homogeneous  system  analytic  in  the  finite  plane  at  a 
sufficient  distance  from  the  negative  real  axis,  with  its  singularities  at  those 
of  dij  (x) ,  b  ( x  +  1 ) ,  and  the  zeros  of  A  ( x ) ,  or  points  congruent  on  the 
left.    It  thus  corresponds  to  the  first  principal  solution  for  /x  >  0. 

These  solutions  will,  as  before,  be  asymptotic  in  the  left  and  right  half 
planes,  respectively,  to  the  series  formally  satisfying  the  equations,  the  first 
of  them  maintaining  its  asymptotic  form  in  the  right  half  plane  at  a  sufficient 
distance  from  the  positive  real  axis.  The  formal  series  have,  however,  a 
somewhat  different  form,  and  now  are 
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gt(x)  =  *r+'l(jr*  +  S^  +  2,  ...,n-l), 


gn)  (x  =  ^(^  +  f-)+  •••)• 


For  if  we  substitute  expressions  of  this  form  in  the  equations,  and  compare 
coefficients  of  like  powers  of  x,  it  is  evident  that  at  each  stage  we  shall  have 
n  equations  in  which  there  are  n  unknown  coefficients  and  coefficients  pre- 
viously determined.  Thus  the  coefficients  can  be  calculated  uniquely.  In 
particular  we  have 

Qi  =  ain  b  (»  =  1,  2,  ••■,7i-l) 

9  =  b. 

A  case  of  interest,  and  one  which  will  be  used  later  in  the  treatment  of  the 
equation  of  the  nth  order,  is  that  where  a#  (x)  =  0,  for  j  =}=  i  +  1,  i  =f=  n.  In 
that  case  the  formal  series  become 

(49)  gi(x)  =  .T"+('">(j7i  +  ^-+  2,  •••,  n-1), 

where 

Qi  —  0>i,  i+l  i+2  '  '  '  nb  (  %  =  1  ,  2  ,   •  •  •  ,  «) 

(50) 

?»=&• 

We  can  also  show  the  existence  of  solutions  for  the  case  y.  <  0  in  a  second 
way  by  referring  it  to  the  case  \x  >  0 .  When  we  remember  the  definition  of 
A  (x) ,  Aji  (x) ,  the  system  (17)  can  be  written 

n 

Qi  {x)  =  £  ai}  (x)     (x  +  1)  +  h  (x  +  1)   ( i  =  i ,  2 ,  •  ■  •  ,  n ) , 

where 


4  ' 


dij  ( x )  =  x~*  I  dj  +      +  •  •  •  ) ,       ai}  =  ( —  1 

&<(*)  =  »v"'t(^  +  ^-+  "•).  (-  Dn+i 

If  now  we  put  x  =  —  (s  +  1 ) ,  gr£  (x)  =  #£  (  —  (z  +  1 )  )  =  (ji  (z  +  1 )  , 
and  ju  =  —  X  (X  >  0 ) ,  this  system  can  be  written 

n 

^(*+l)=Eo</(a)&(*)  +  6i(*+l)   (i-l,  2,  ....  n), 

7=1 
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where,  when  \z  \  >  R, 


oi3  =  (-  iy+i^f, 


Ani  b 
A 


Since  X  is  positive  the  previous  existence  theorems  can  be  applied.  The  sub- 
stitution x  —  —  ( z  +  1 )  shows  that  the  roles  of  right  and  left  half  planes 
have  been  interchanged. 

We  can  summarize  the  above  results  as  follows: 

When  p  <  0 ,  first  and  second  'principal  solutions  exist  and  have  their  singu- 
larities at  the  same  points  as  before.  The  first  solution  is  asymptotic  to  the  formal 
series  in  the  right  half  plane;  the  second  solution  is  asymptotic  to  the  same  series 
throughout  the  left  half  plane,  and  also  in  the  right  half  plane  at  a  sufficient  distance 
from  the  real  axis. 


It  still  remains  to  consider  the  question  of  solutions  for  the  intermediate 
case  p  =  0 .    It  will  appear  that  at  times  this  case  is  exceptional. 

We  shall  first  examine  when  the  solutions,  as  found  on  the  hypothesis  that 
p  >  0,  are  valid  when  we  put  p  =  0 .  Consider  the  equation  (15),  and  let 
us  see  under  what  conditions  the  series  will  still  converge.  The  typical  term 
can  now  by  (16)  be  written 


If  |  pj  [  >  1 ,  we  shall  have  a  series  in  which  the  terms  decrease  in  practi- 
cally geometrical  fashion,  so  that  the  series  converges  uniformly.  Since 
\Pn  |<  I  Pi  \  {i  <  n) ,  the  series  for  coi  (x) ,  •  •  •  ,  co„  (x)  will  all  converge  if 
|  pn  |  >  1 .  The  first  principal  solution  as  found  for  p  >  0  therefore  still 
exists  for  p  =  0 ,  if  |  pn  \  >  1 .  Under  this  same  condition  the  second  principal 
solution  for  ju  >  0 ,  given  by  the  contour  integrals,  is  also  valid.  For  in  this 
case  (39)  can  be  written 


§  10.    The  case  p  =  0 . 


(x  +  m  )"-'' 


[  df,  x  -f  to  ] . 


and  this  expression  decreases  along  A  <x>  and  B  <x>  in  such  a  way  that  the 
integral  in  (36)  converges. 
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We  can  then  say:  If  every  root  of  the  characteristic  equation  is  greater  than 
unity  in  absolute  value,  we  can  consider  the  case  p  =  0  as  a  limiting  case  for  p  >  0 . 

In  a  similar  way  we  can  show  that:  When  all  the  roots  of  the  characteristic 
equation  are  less  than  unity  in  absolute  value,  the  case  p  =  0  can  be  considered 
as  a  limiting  case  for  p  <  0 . 

We  still  have  to  determine  whether  we  can  make  the  sum  formulae  valid 
if  (1)  some  of  the  roots  of  the  characteristic  equation  are  greater  in  absolute 
value  than  unity  and  others  less,  or  (2)  some  of  the  roots  of  the  characteristic 
equation  are  equal  to  unity  in  absolute  value. 

Consider  the  first  case.  Let  pi ,  •  •  •  ,  p&  be  greater  than  unity  in  absolute 
value,  and  pk+\ ,  ■  ■  ■ ,  pn  he  less.  Then,  if  we  take  as  solutions  of  the  homo- 
geneous equations  the  first  principal  solutions,  we  can  use  for  the  sum  formulae 
for  coi  (  x  )  ,  •  •  •  ,  ojfc  ( x )  the  direct  sum  to  the  right.  We  will  next  show  how  to 
express  Uk+i  ( x ) ,  •  •  • ,  un  ( x )  as  contour  integrals.  As  the  work  is  similar 
for  all  of  them,  we  shall  consider  the  equation 

,5!)  ABW(,)-(-l)^t-('+^f+l)!,+  1)- 

If  x  is  in  the  right  half  plane  and  is  sufficiently  large,  we  know  by  (16) 
that  the  right-hand  side  of  (51)  can  be  written 

 ^  [ dk+i;  x+  1  J. 

PiH-i 

In  §  1  it  was  stated  that  there  is  a  set  of  intermediate  solutions  of  (2)  that  are 
asymptotic  to  the  formal  series  in  the  right  half  plane  with  regard  to  x,  and 
with  regard  to  v  in  the  left  half  plane.  Suppose  these  were  the  solutions  of 
the  homogeneous  system  that  we  were  using;  we  could  then  write 

(51')         AwiH-i  (*)  =  (-  i)n+k+l(x+xly  Wi;  *+  1], 

where  [du+i,  x  +  1  ]  is  bounded  for  x  large  in  the  right  half  plane  and  for  v 
large  in  the  left  half  plane.  Since  |  pk+i  |  <  1 ,  it  is  evident  that  for  v  large 
enough  the  above  value  of  Aco^+i  (x)  decreases,  when  x  goes  to  infinity  in  a 
direction  parallel  to  the  negative  real  axis,  approximately  as  an  exponential. 
The  equation  could  then  be  solved  by  means  of  a  contour  integral  sum  to 
the  left.  Return  now  to  the  consideration  of  equation  (51).  The  first  prin- 
cipal solutions  can  be  expressed  in  terms  of  these  intermediate  solutions  by 
equations  similar  to  (8),  and  reference  to  (37)  and  (38)  shows  that  the  right 
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hand  side  of  (51)  would  be  replaced  by  a  group  of  terms  such  as  the  one  on 
the  right  in  (51')  in  which  every  p  is  less  in  absolute  value  than  unity.  From 
this  we  see  that  o)k+i  ( x )  can  be  expressed  as  a  contour  integral  about  a  con- 
tour to  the  left,  and  similarly  coa-+2  ( x  )  ,  ■  •  ■  ,  co„  ( x ) . 

If,  therefore,  we  use  the  first  principal  solutions  of  the  homogeneous  system, 
and  divide  the  resulting  sum  formulae  into  two  groups,  according  as  they  correspond 
to  roots  of  the  characteristic  equation  greater  in  absolute  value  than  unity,  or  less, 
we  can  express  the  first  group  by  a  direct  sum  to  the  right,  and  the  second  by  a 
contour  integral  to  the  left.  Thus  we  have  the  first  principal  solution  of  the 
non-homogeneous  system. 

In  a  similar  way  we  can  show  that  if  we  use  the  second  principal  solutions 
of  the  homogeneous  system,  and  divide  the  sum  formulae  as  before,  we  can 
solve  the  first  group  by  a  contour  integral  to  the  right,  and  the  second  by  a  direct 
sum  to  the  left.    This  gives  a  second  principal  solution. 

The  only  remaining  case  is  that  in  which  some  of  the  roots  of  the  characteristic 
equation  are  equal  to  unity  in  absolute  value.  This  case  has  not  been  com- 
pletely investigated  for  the  homogeneous  system.  Here,  in  certain  cases,  at 
least,  a  modified  sum  may  give  a  convergent  series.  An  example  of  this  is 
the  equation  occurring  in  the  Weierstrass  development  of  the  gamma  function, 
namely  the  equation 


is  also  a  solution,  and  converges. 

When  ju  =  0  the  asymptotic  form  of  the  solution  requires  a  slightly  different 
treatment.  Take  the  case  where  all  the  roots  of  the  characteristic  equation 
are  greater  than  unity  in  absolute  value,  and  consider  the  asymptotic  form  of 
the  first  principal  solution. 

We  need  only  find  the  asymptotic  form  of  co,  ( x )  hij  ( x ) ,  where 


g  (x+  1)  -  g  (x)  = 


1 


x 


The  direct  sum  to  the  right  is 


which  diverges,  while  the  modified  sum 


m—l 


(x+m)"-rj 
Pj+m 


[d,-;  x  +  m}. 
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Take  the  first  k  terms  of  wj  (x) ,  where  J  |  x  \  —  1  <  k  <  J  |  x  |;  we  can  write 
them  in  the  form 


[  dj;  x  +  m  ] . 


In  all  these  terms  the  expression  (1  +  m  j  x  )v~ri  can  be  developed  in  a  series 
in  1  /  x ,  since  |  m  /  x  |  <  1 .  There  is  a  marked  difference,  however,  between 
these  developments  and  those  for  the  case  /x  >  0  as  given  in  §  4.  In  the 
present  case  each  of  the  k  terms  of  the  summation  when  expanded  will  contain 
a  constant  term,  while  in  the  corresponding  expression  (20)  in  the  former  case, 
the  first  term  alone  when  expanded  contained  a  constant  term,  the  others 
beginning  with  aT" ,  •  •  •  ,  aT(fc-1)M .    The  final  constant  term  in  the  present 

h 

instance  will  evidently  be  dj  zl  pjm .    When  k  increases  this  approaches 

m=l 

dj/  (pj  —  1 ) ,  since  |  pj  |  >  1 .  The  first  k  terms  of  co;-  ( x ) ,  when  multiplied 
by  ha  (x) ,  can  therefore  be  written  a;"[Xjft);  x],  where  \j°  approaches 
dj  S{j  /  ( pj  —  1 )  as  k  increases. 

The  terms  of  coy  ( x )  after  the  kth  behave  ultimately  like  a  geometrical  pro- 
gression, since  the  ratio  of  successive  terms  obviously  approaches  1  /  pj  as 
the  number  of  the  terms  becomes  large.  It  follows  then  that  the  remainder 
after  k  terms  in  the  series  for  co,-  ( x ) ,  is,  provided  k  is  large  enough,  of  essenti- 
ally the  order  of  the  (&+  1  )th  term,  which  is  M  (x+k+1  y-r'/p*+k+1 ,  where 
M  is  bounded.    The  product  of  this  and  ha  (x)  can  be  written 


where  |  M  |  is  less  than  some  fixed  constant.  Since  k  increases  in  a  constant 
ratio  to  x ,  this  becomes  smaller  than  any  power  of  1  /  x ,  so  that  we  have  ulti- 
mately aij  ( x )  hij  (x)  =  xv  [~Kjk);  x]  where  Xj4)  approaches  djSal  (pj  —  1)  for 
k  large. 

In  this  way  we  can  determine  the  asymptotic  form  of  the  solution,  and  we 
are  led  to  try  to  find  formal  solutions  of  the  form 

gi(x)  =  »V(^  +  ^-+  2,  •••  ,  n). 

That  the  coefficients  in  these  expressions  can  be  uniquely  determined  so  as 
to  satisfy  the  equations  can  be  shown  easily  by  substitution.  In  particular 
we  find 

bA  ■ 

(52)  gi  =  (-  l)***H-=?  (<  =  1,  2,  ...  ,  n), 

A 
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where  A  is  the  determinant  |  a#  —  5#  \ ,  ( 5#  =  0 ,  if  i  4=  j ,  8u  —  1 ) ,  and  Ai 
is  the  minor  of  A  with  reference  to  the  element  in  the  last  row  and  iih  column. 
Also  we  find  again  that  the  first  and  second  principal  solutions  are  asymptotic, 
in  the  right  and  left  half  planes  respectively,  to  the  series  formally  satisfying  the 
equations. 

§  11.    The  equation  of  the  nth  order. 

We  shall  consider  a  single  non-homogeneous  equation  of  the  nth  order, 
which  we  write  in  the  form 

(53)    g  (x  +  n)  +  a,  (x)  g  (x  +  n  -  1)  +  •  •  •  +  an  (x)  g  (x)  =  b  (x) . 

The  functions  di  ( x ) ,  ■  •  •  ,  an  (x) ,  b  (x)  are  assumed  to  be  rational,  so  that 
we  have 

,(D 


/        0a>  \ 

ai(x)  =  sfri  ai  +  ^rH  ), 

(     bw  \ 

b(x)  =  xv\b  +  —  +'■'), 


(for  |  x  |  >  R  )  . 


In  place  of  the  single  equation  (53),  consider  the  system 

<7i  (2C+  1)  =  x»  g2  (x), 
g2  (a-+  1)  =  .rM  g3  (x) , 


gfn_!  (cc+  1)  =  x*  gn  (x), 


(54)     gn(x+  1)  =  —  ?T7  ;_,w  TV  7  „_  iW^C*) 


„v— (n— 


+ 


/       bw  \ 


If  from  these  equations  we  eliminate  #2  ( £ ) ,  •  •  •  ,  gn  ( .r ) ,  we  shall  find  that 
gx{x)  satisfies  (53),  so  that  the  system  is  equivalent  to  the  original  equation  of 
the  nth.  order.    When  we  expand  the  quantities  ( 1  -f-  jjxY  ,  which  occur  in 


■ 
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the  last  equation  above  we  find  that  that  equation  can  be  written 

gn  (x +  1 )  =  Si  (x )  gfi  (cc)  +  S2  {x )  g2  ( a)  +  •  •  •  +  o„  ( .r )  gr„(  .t )  +  &  ( a-  +  1 ) , 

where 


Cli  (x)  =  —  X*  (  Cfn_£+1  +  ~  + 


(»+-+■•■•), 


so  that  the  system  is  of  precisely  the  form  (1).  The  determinantal  charac- 
teristic equation  reduces  to 

(55)  pn  +  oi  p"-1  +  h  «„-i  p  +  a„  =  0 , 

which  can  therefore  be  found  directly  from  (53). 

We  can  now  apply  to  the  system  the  existence  theorems  already  estab- 
lished, and  we  thus  have  solutions  gi  ( x )  of  (53)  in  all  cases  except  when  /x  is 
zero  and  some  of  the  roots  of  the  characteristic  equation  are  equal  to  unity  in 
absolute  value. 

The  asymptotic  form  of  the  solutions  of  (54)  can  be  determined  by  applying 
the  general  results  for  a  system  (for  /x  <  0,  apply  (49)).  In  each  case  it  is 
found  that  the  solution  gx  ( x )  is  asymptotic  to  the  series  formally  satisfy- 
ing the  equation  (53).    This  formal  solution  is  now 

g  (x)  =  x"'n>i  (^g  +  ^7  +  •  •  ( for  M  >  0  )  , 

g  (x)  =  xV{j3  +  9—+  ■•')  (for/x^O), 


where 


b 

g  =  —  ( for  ai  >  0  )  , 


g  =  b  (  f  or  ix  <  0  ) 

a  =  —  ■  —  :   ( for  £1  =  0). 

1  +  ai  +  a2  +  •  ■  •  +  an 

The  constant  term  can  in  all  cases  be  found  either  from  the  known  form  of 
the  constant  term  of  the  solution  of  the  corresponding  system  (for  /x  <  0, 
apply  (50);  for  /x  =  0,  use  (52)),  or  by  direct  substitution  in  equation  (53). 
We  can  now  state: 
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Theorem  IV.  If  fx  <  0 ,  and  all  the  roots  of  the  characteristic  equation  (55)  are 
different  from  zero,  the  equation  (53)  has  two  principal  solutions  gi(x)  and  g{  (x) . 
The  solution  gi  (x)  is  analytic  in  the  finite  plane  save  for  possible  poles  at  the 
poles  of  ai  ( x ) ,  •  •  •  ,  a„_i  ( x ) ,  b  ( x ) ,  at  the  zeros  of  a„  ( x ) ,  or  at  points  con- 
gruent on  the  left;  while  g{  (x)  is  analytic  save  for  poles  at  the  poles  of  at  (x) , 
•  •  •  ,  an  ( x ) ,  6  ( x )  or  points  congruent  on  the  right.  If  [i  =  0 ,  same  is 
£r«e  provided  in  addition  that  none  of  the  roots  of  the  characteristic  equation  are 
equal  to  unity  in  absolute  value.  There  exists  a  series  g  ( x )  formally  satisfying 
(53),  and  in  all  cases  in  which  ive  have  shown  the  existence  of  solutions  we  have 
<7i  (x)  ~  g  (x)  in  the  right  half  plane,  while  g\  (x)  ~  g  (x)  in  the  left  half 
plane. 

Princeton  University, 
May,  1912. 


